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We utilize a discrete (sequential) measurement protocol to investigate quantum process tomog¬ 
raphy of a single two-level quantum system, with an unknown initial state, undergoing Rabi oscil¬ 
lations. The ignorance of the dynamical parameters is encoded into a continuous-variable classical 
system which is coupled to the two-level quantum system via a generalized Hamiltonian. This com¬ 
bined estimate of the quantum state and dynamical parameters is updated by using the information 
obtained from sequential measurements on the quantum system and, after a sufficient waiting pe¬ 
riod, faithful state monitoring and parameter determination is obtained. Numerical evidence is used 
to demonstrate the convergence of the state estimate to the true state of the hybrid system. 

PACS numbers: 03.65.Ta, 03.65.Wj, 06.20.-f 


I. INTRODUCTION 

Recent years have witnessed remarkable progress in the 
theoretical study m and experimental realization P- 
E] of the manipulation of single quantum systems in 
the framework of quantum information processing and 
communication. The ability to monitor and control sin¬ 
gle quantum systems [IDHI3] is essential for the devel¬ 
opment of technology based on these quantum building 
blocks, hence methods for quantum-state monitoring and 
dynamical parameter estimation are of paramount im¬ 
portance. 

A novel method for tracking the evolution of a driven, 
isolated two-level system in real-time by means of a se¬ 
quence of measurements with minimal disturbance, was 
devised by Audretsch et al. [HHis] and is related to 
the concept of continuous measurement Pi] and state 
estimation m- The state estimate (or state guess) 
and the real state of the quantum system are evolved 
with the same propagator given by the Hamiltonian and 
measurement-dependent dynamics. It was argued that, 
for known dynamics but an unknown initial state, the 
estimated state and real state eventually converge m- 
Remarkably, numerical simulations show that the con¬ 
vergence for two-level systems undergoing Rabi oscilla¬ 
tions occurs within only a few Rabi cycles [H] . The the¬ 
ory also applies to many-level systems (such as a Bose- 
Einstein condensate in a double-well potential m) and 
systems with inhnite-dimensional Hilbert spaces (for ex¬ 
ample, particles in various potentials m)- Similar meth¬ 
ods have been used to experimentally control the number 
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of photons in a cavity m, the frequency of Rabi oscil¬ 
lations in a superconducting qubit m, and for the state 
tomography of an ensemble of cold cesium atoms p. 

However, these methods fail when the dynamics of the 
system are not precisely known - for example, there may 
be some uncertainty in the Rabi frequency for an atom 
driven by an external laser field. Quantum process to¬ 
mography, which is the task of identifying the dynami¬ 
cal quantum process, now becomes necessary. The stan¬ 
dard method is extremely resource intensive since the 
dynamical parameters are inferred through the measure¬ 
ment of an informationally complete set of observables, 
at different times, on a large ensemble of identically pre¬ 
pared quantum systems [20j . As an efficient alternative, 
the techniques of direct characterization of quantum dy¬ 
namics (DCQD) 1^ and selective and efficient quantum 
process tomography (SEQPT) were devised. DCQD 
substantially reduces the number of resources required 
for process tomography by utilizing quantum correlations 
between the probe system and an ancillary qubit. SE¬ 
QPT, on the other hand, achieves this objective by map¬ 
ping the estimation of the dynamical parameters (speci¬ 
fied a process matrix) to the average fidelity of the quan¬ 
tum channel which can, in turn, be estimated by ran¬ 
domly sampling over a particular set of states called a 
“2”-design. However, these methods use projective mea¬ 
surements and therefore destroy the coherent dynamics. 

In contrast, the sequential measurement protocol in¬ 
troduced here allows one the possibility of sensing the 
dynamical parameters in real time by measuring a sin¬ 
gle quantum system as its state evolves. The incomplete 
knowledge of the parameters is encoded into the state of 
a hypothetical classical system, which is coupled to the 
quantum system via a generalized Hamiltonian. Both 
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systems are combined to form a hybrid quantum-classical 
system. By means of this construction, the measure¬ 
ment record obtained from sequential measurements on 
the quantum system can be used to update the estimated 
hybrid state and thus the knowledge of the dynamical pa¬ 
rameters. Essentially, we demonstrate that the problem 
of parameter estimation can be reduced to that of state 
estimation on a higher-dimensional state space. Numer¬ 
ical simulations demonstrate that the convergence of the 
estimated hybrid state to the real state is achievable in 
real time. It is quite remarkable that the quantum state 
need not be known in order to determine the dynamical 
parameters. 

Attempts to derive an estimation scheme for the dy¬ 
namical parameters via continuous measurements were 
made by Negretti and Mplmer , Ralph et al. [21] , and 
very recently by Six et al. [53]. However, they used a 
different estimation method and derived separate update 
equations for the quantum system and probability distri¬ 
bution of the parameters. 

The remainder of the paper is structured as follows. 
Section |Il| briefly introduces the notion of unsharp mea¬ 
surements as an essential tool for quantum-state estima¬ 
tion and monitoring. Section III describes how to extend 
the state-estimation method for the determination of un¬ 
known parameters within the framework of the hybrid 
quantum-classical theory. Numerical results are analyzed 
in Sec. US while Sec. 0 contains a concise discussion of 
the estimation procedure. Section |Vl| concludes. 


II. UNSHARP MEASUREMENTS AND STATE 
ESTIMATION 

We are interested in the monitoring, or continuous ob¬ 
servation, of the evolution of a single quantum system in 
real time. Consider a two-level quantum system under¬ 
going Rabi oscillations due to the influence of the Hamil¬ 
tonian 


( 1 ) 

where is the Rabi frequency and ux is the Pauli ma¬ 
trix that generates rotations about the X axis. This sys¬ 
tem could be, for example, two levels in the hyperfine 
ground-state manifold of a trapped ion driven by a mi¬ 
crowave field. The normalized state is given by 


also be used when we have only a single quantum sys¬ 
tem, but requires many repetitions of the experiment for 
each delay time t. As an alternative, we employ an es¬ 
timation method that uses a sequence of a special kind 
of positive-operator valued measure (POVM) measure¬ 
ments; so-called unsharp measurements | 26 j . in order to 
gain information about the state of the system with min¬ 
imal disturbance. 

A POVM measurement with outcome n, on a system 
in the state I'i/'), will result in the postmeasurement state 

MnW (3) 

where Mn is the Kraus operator corresponding to the 

measurement result n. In order to monitor the dynamics 
of the oscillating two-level system we perform unsharp 
measurements of the observable. The Kraus operators 
are thus given by 

^0 = \/l^^|0)(0|-k v^|l)(l| (4) 

Mi = y^\0){0\ +^/l-po\l){l\, ( 5 ) 

related via MqMq + M\Mi = I and 0 < po < 0 - 5 - The 
strength of the discrete measurement sequence depends 
on the strength of the individual measurements Ap = 
Po — (1 — po) = 2po — 1, as well as the frequency of the 
measurements 1/r. It is best quantified by the level- 
resolution rate 7m = 1/7V„ with = T/(Ap)^ [7|. 

The estimation and monitoring method of Diosi et 
al. m now proceeds as follows. A sequence of unsharp 
measurements, with a frequency l/r, is applied to the 
two-level system as it evolves in time according to the 
known Hamiltonian. The duration of the measurement 
is assumed to be much smaller than all other dynami¬ 
cal timescales (impulsive measurement approximation), 
therefore the state change due to measurement can still 
be described by Eq. (|^. After N measurements at time 
t = Nt, the system is in the state 

|V> (Nt)) = Mr^^U (Hfl, (Hk, t) ... 

X (6) 

up to an appropriate normalization constant, where 



U {^r,t) = exp 


~H {VLr)t 


(7) 


\ilj{t))=co{t)\0}+ci{t)\l), (2) 

where |cop-l-|cip = 1 and |0) (|1)) represents the ground 
(excited) state of the system. The usual way to determine 
the dynamics (or Rabi oscillations) of |ci(f)p involves 
the preparation of the initial state of the quantum sys¬ 
tem, followed by time-evolution of the coherent dynamics 
for some delay time t, ending with a projective measure¬ 
ment. This procedure is repeated for different delay times 
with an ensemble of quantum systems. The method can 


The same sequence of operators is then applied to a pure 
estimated (or guessed) state lipe), which is orthogonal to 
the real state in the worst case. It has been argued in 
Ref. [17] that the effect of a sequence of measurements 
and measurement-independent unitary evolutions on the 
state of the system can be approximated in the limit 
of continuous measurements. In this limit the evolution 
is described by coupled Ito stochastic master equations 
for the true state of the system \tj}) of the system, the 
estimated state \^pe) and the measurement record. The 
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analytical methods of stochastic Ito calculus were used to 
show the convergence of the estimated state to the real 
state under continuous evolution and measurement m- 
Moreover, the convergence of the estimated state has 
been tested by numerical simulations for various systems 

[nmmis]. 


III. PROCESS TOMOGRAPHY 

The aforementioned method of state estimation and 
monitoring only works accurately if all parameters of 
the Hamiltonian are known precisely. Although it allows 
monitoring with high fidelity in the presence of small con¬ 
tinuous noise m and infrequent scattering events m, 
a lack of knowledge of any of the dynamical parameters 
may result in completely unfaithful outcomes. 

In the following, we consider monitoring the state 
of a system with a single unknown dynamical param¬ 
eter; for example, the Rabi frequency. In general, we 
know the Rabi frequency to be in some finite range, i.e., 
G [wminjWmax]- For simplicity, we assume that there 
are N discrete frequencies, uji, within this range and the 
real frequency H/j is one of them. Initially, we assign 
an equal probability P{uji) (i.e., an unbiased probability 
distribution) to the frequencies. This corresponds to hav¬ 
ing no knowledge about which of the frequencies is the 
correct Rabi frequency. The first step of state estima¬ 
tion involves propagating the estimated state according 
to the Hamiltonian dynamics for a time period r. How¬ 
ever, under unitary time evolution [Eq. 0] with an un¬ 
known frequency, the initially pure estimated state |i/'e) 
naturally evolves into a mixed state 

Mt) = '^ P{^i)UiuJi,T)\ll;e){lpe\U\uJi,T). ( 8 ) 

i 

We will now demonstrate how the frequencies LOi [spec¬ 
ified in [/(a;i,r)] can be incorporated into an effective 
higher-dimensional state so that the technique of state 
estimation can still be applied. 

We aim to formulate our estimation method for the 
quantum state and Rabi frequency in terms of the hy¬ 
brid quantum-classical theory |27II29j . For introductory 
purposes, we first describe the estimation technique by 
means of a formalism which allows us to represent the 
unknown dynamical parameter as the state of an addi¬ 
tional quantum system [23l[3D]. In particular, we artifi¬ 
cially construct a diagonal density matrix (i.e., a mixed 
quantum state) to represent the probability distribution 
of the frequencies. We refer to this representation as the 
quantum-quantum formalism to distinguish it from the 
hybrid quantum-classical formalism. The latter reflects 
better the classical nature of the probability distribution 
as the mixed state of a classical system. Although we 
consider only a single unknown parameter, the result can 
easily be generalized to full process tomography, as shown 
in Sec. ITVl 


A. Quantum-Quantum formalism 

The initial probability distribution of the frequencies, 
at time t = 0, can be written as the state of a hypothetical 
quantum system as follows 

Pclass(O) — ^ ^ -F(cl^^) jcU^) (cU^ I, (9) 

i 

where P{^i) = 1 ^ind no off-diagonal elements (co¬ 
herences) are allowed for pciass in this formulation [23j . 
The combined estimate of the frequency and the two- 
level state (at time t = 0) can then be expressed as the 
tensor product of the two corresponding states 

/ 3 ''( 0 ) 0 ( 10 ) 

i 

Such a density operator acts on the Hilbert space 
where S^uj is the Hilbert space spanned by the orthonor¬ 
mal vectors \uji) corresponding to all possible dynamical 
parameters Wi and is the Hilbert space for two-level 
quantum systems. We can also upgrade the Hamiltonian 
[in Eq. Q] in a similar fashion and redefine it as 

n = ^\uJi){uJi\iSi (II) 

i 

The estimated state p® [Eq. ®] evolves under this 
Hamiltonian as 

='^P{uJi)\uJi){uJi\ U{uJi,T)\lpe){fpe\U\uJi,T). 

(12) 

If we take the partial trace with respect to the first 
system, then we obtain the quantum state defined by 
Eq. Hence, this construction allows us to shift the 
ambiguity in the dynamical parameters from the Hamil¬ 
tonian to the composite state. The task is now to esti¬ 
mate the Rabi frequency and the state of the two-level 
quantum system using the composite state. 

The estimation experiment proceeds as previously de¬ 
scribed in Section |llj The quantum system evolves un¬ 
der the Hamiltonian for a time r after which a 

single unsharp measurement is performed on it, yielding 
measurement result n. The estimated composite state si¬ 
multaneously evolves under H for the time r and, based 
on the measurement result, the augmented measurement 
operator 

7W„=I0M„. (13) 

is applied on it. An unsharp measurement on the quan¬ 
tum system yields information not only about the true 
state of the system but also about the Rabi frequency, 
and after an evolution for time r and a single measure- 
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merit, the estimated composite state is updated as follows 
p®(r) p®'(t) 

^ ' i 

^ MnU{uj,,TMe){i’e\U\i0,,r)Ml (14) 

= -^-r'^P{uJi)P{n\uJi)\uJi){uJi\ (g)/3e(Wj,T), 
r(n) 

(15) 


where 


Pe(w*,r) = 


P(n|wi) 


M„ (/(wi, r) I V'e) (V’e IP’'' (wi, r) , 


(16) 


P{n\uJi) = Tr [M„;7(w„r)|V',)(r^,|;7t(cc„r)Mt] , (17) 
P{n) = Y. P{ijJi)P{n\uJi). (18) 


P{n\uJi) is the probability of measuring result n for the es¬ 
timated state under the condition that it evolved through 
the Hamiltonian H{uii) = huJiCTxl‘2. up to a time r, when 
the measurement was performed on the system. After 
each measurement, the observer’s knowledge about the 
Rabi frequency is thus updated as follows, 

P{uji) P{uji\n) = -^—P{n\u}i)P{uJi). (19) 

P(n) 

This resembles the update of probabilities according to 
Bayes’ Law [5T| . 

As previously mentioned, after a large number N of up¬ 
dates, the state-estimation method leads to convergence 
of the estimated state to the true state of the system. 
Since we have specified that the true Rabi frequency 
is one of the frequencies uii, the probability density will 
eventually approach a Kronecker-Delta function indicat¬ 
ing that the correct frequency has been determined, i.e., 

p'^{Nt) g |V’e(^'r))(lAe(^T)| 

i 

= g \'lpeiNT)){'lpeiNT)\. (20) 


The estimation fidelity, which measures the overlap be¬ 
tween the real and estimated quantum states, is exactly 
unity after this time, demonstrating perfect state moni¬ 
toring of the single quantum system in real time, as well. 


B. Hybrid formalism 

We now consider estimation of the Rabi frequency and 
the state of a two-level quantum system within the the¬ 
ory of hybrid quantum-classical systems by translating 
the method described in the previous section. In the 
hybrid formalism, the probability distribution of the fre¬ 
quencies, P{uji), can be viewed as the statistically mixed 


state of a hypothetical classical system while the estimate 
of the quantum state is treated quantum mechanically. 
For the estimation procedure, we construct a hybrid sys¬ 
tem, whose hybrid state (at time t = 0) is given by 

p^{u}i,0) = P{uJi)\ll^e){lpe\- ( 21 ) 

This state is positive semidefinite 

P'’(wi,T)>0, i = l,2, ...,A^ (22) 

and normalized 


Tr 




= 1 . 


(23) 


The real state of the quantum system (at time t = 0) 
can be written as 


p(0) = (24) 

which evolves under unitary dynamics as 

p{T) = u{nR,T)mu\nR,T). (25) 

We prescribe the following evolution for the hybrid state, 

= U{uJi,T)p={uj„0)U\uJi,T). (26) 

The uncorrelated (product) structure of the hybrid state 
[Eq. pT] )] is immediately lost after this operation. In the 
hybrid formalism, the sum over the frequencies reveals 
the reduced state of the quantum subsystem, i.e., 

Pe('r) = (27) 

i 

which is precisely what is specified in Eq. (|^. On the 
other hand, the trace over the hybrid state gives the re¬ 
duced state of the classical system 

P(a;,) =Tr[p-(cc„T)] . (28) 

Due to a measurement with outcome n, the real quan¬ 

tum state changes like 

p{t) hO M„p(t)M,( (29) 


and we update the estimated hybrid state as follows: 

p^{uJi,T) (30) 

The joint probability for the frequency uji and the mea¬ 
surement result n can be obtained from the hybrid state 
as 


P{oJi,n) 


Tr 


MnP^{uJi,T)Ml 


(31) 


The updated probability distribution of the frequencies 
thus takes the equivalent form: 


P(a;i|n) 


P{^^i,n) 

Ei^(wi,n) 


1 

P(n) 


P{n\uJi)P{uji). 


(32) 


This result coincides completely with the quantum- 
quantum result, cp. Eq. (19). 
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IV. NUMERICAL SIMULATIONS 
A. Frequency estimation 

We can now test the performance of our state and pa¬ 
rameter estimation method via numerical simulations. 
We want to determine the frequency of the Rabi oscil¬ 
lations for a two-level ion due to an external driving 
field as well as monitor the state of the ion. Let us as¬ 
sume that the Rabi frequency is known to lie in the 
range [0.95r2o;l-05flo] where Oq is an experimentally de¬ 
termined value. This corresponds to a relative frequency 
error of 5% and we wish to track the Rabi frequency with 
a relative error of 1%. We therefore use a discrete grid 
of values in the specified range where the spacing of the 
values on the grid is Oq/IOO, which is the required ac¬ 
curacy for the frequency. An unsharp measurement of a 
single observable, such as in our case, on the quan¬ 
tum system is sufficient to determine the frequency. For 
our sequence of measurements, we select Ap = 0.2 as 
the strength of the individual measurements and a mea¬ 
surement period of r « Tr/10, where Tr = 27r/nR is the 
Rabi period. This selection ensures that the Rabi oscilla¬ 
tions are only weakly disturbed due to the measurement 
sequence since the measurement strength « IIr/(Stt), 
is smaller than the Rabi frequency. This is important for 
monitoring the state of the ion, since a strong measure¬ 
ment would immediately project the ion into either the 
ground or excited state, i.e., it would freeze the dynamics 
(similar to the quantum Zeno effect [5^1. 

Since we assume no knowledge of the actual Rabi fre¬ 
quency within the specified range we choose an initially 
flat (or unbiased) probability distribution of the frequen¬ 
cies while the actual probability distribution is taken to 
be a Kronecker-Delta function at the correct frequency. 
The initial quantum state estimate is chosen to be orthog¬ 
onal to the real state, which is the worst case scenario. 
We perform 5000 measurements (or approximately 500 
Rabi cycles) averaged over 1000 runs. This computation 
takes only 550 seconds on a desktop computer with a dual 
core processor. In addition, we provide, in Section |Vj a 
method to approximately halve the computational time 
required for the determination of the frequency. Fur¬ 
ther optimization techniques will be discussed in a future 
work. 

In Fig.j^we plot the classical estimation fidelity for the 
probability distribution of the frequencies as a function of 
the number of measurements for a single run (dashed blue 
line) and averaged over 1000 runs (red line). The estima¬ 
tion fidelity is calculated by using the well-known formula 
fo'' probability distributions pi and qi. The 
graph asymptotically tends to unity and at this point we 
are able to determine the frequency within the specified 
accuracy. An estimation fidelity for the frequency with 
arbitrary precision can be obtained by decreasing the grid 
point distance. If the relative error with which we track 
the frequency is smaller than 0.1% then we would also 
achieve perfect state monitoring once the frequency has 



FIG. 1. (Color online) Estimation of the Rabi frequency. The 
estimation fidelity for the probability distribution of freqnen- 
cies as a function of the nnmber of measurements for a single 
run (dashed blue line) and averaged over 1000 runs (red line). 


been determined m- 


B. Full process tomography 


Consider the case of a spin-^ particle precessing due 
to a magnetic field. The Hamiltonian is given by 


Hi^L) = • CT 


with 


= J2n,c 


(33) 


where Hr is the Larmor frequency, Gi are the Pauli ma¬ 
trices for the i = x,y,z axes and the normalized vector 
n represents the direction of the magnetic field. If the 
magnitude and direction of the magnetic field are not 
precisely known, then there will also be an uncertainty 
in the Larmor frequency and the axis of rotation of the 
particle. In order to fully determine the rotation axis, it 
is not possible to gain sufficient information from mea¬ 
surements of an observable along a single axis. We must, 
therefore, utilize informationally-complete (IC) unsharp 
measurements of noncommuting orthogonal spin observ¬ 
ables. We can construct the required measurement oper¬ 
ators from the projectors which measure the spin around 
the X, y and z axes, Pj_ = | (I ± &{) with i = x,y, z. The 
measurement operators are thus given by 

^ (yr^p; + V^Pl) (34) 

Ml = ^ (35) 

related via Mq^Mq + Mf Mq -I- Mf + 

= I 

For the simulation we consider the situation where the 
Larmor frequency has a relative error of 5% and the axis 
of rotation (specified by a Bloch vector with parameters 
9 and c^) is completely unknown. We select 10 points in 
the range [0.95Ho,1.05Ho] for experimentally determined 
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Number of Measurements 


FIG. 2. (Color online) Estimation of the Larmor frequency 
and axis of rotation for a spin-| particle. The estimation 
fidelity for the probability distribution of the dynamical pa¬ 
rameters as a function of the number of measurements for a 
single run (dashed blue line) and averaged over 1000 runs (red 
line). 


Hq 10 points each for the parameters 9 and (p, in the 
ranges [0,7r] and [0,27r], respectively. We once again se¬ 
lect the parameter Ap = 0.2 for the individual measure¬ 
ments and a measurement period of r « T^/IO where 
Ti^ = 27r/ni is the Larmor period. The initial state esti¬ 
mate of the quantum system is taken to be orthogonal to 
the real state. We perform 30000 measurements on the 
quantum system and update the hybrid estimate accord¬ 
ingly. Figure shows the convergence of the estimated 
probability distribution to the actual probability distri¬ 
bution. We can observe that the fidelity tends, asymp¬ 
totically, to unity. 


(leading to the Quantum Zeno effect for continuous pro¬ 
jection measurements) and the original Rabi frequency 
cannot be detected. On the other hand, if the measure¬ 
ments are very weak, i.e., Tm ^ Tfl, the Rabi oscilla¬ 
tions are not modified but there is little information gain 
per measurement and the convergence of the probability 
distribution of frequencies is slow. In practice, there is 
an optimal measurement strength that allows one to re¬ 
solve the dynamical parameters with a certain accuracy 
in a minimum time. For example, here we used mea¬ 
surements with medium level-resolution time « Tr 
to resolve the Rabi frequency and the direction of the ro¬ 
tation axis with moderate accuracy. Optimal sequential 
measurements are the subject of our current and future 
studies, which will be presented elsewhere. 

One particular advantage of our method is that the 
estimate or state guess can be a pure state (represented 
by a state vector) instead of a density matrix. We can 
now utilize this fact to substantially decrease the com¬ 
putational time required to determine the dynamical pa¬ 


rameters. Hence, in the hybrid formalism, Section IIIB 


we replace the density matrix for the estimate state with 
a pure state, i.e.. 


p'=(a;i) ^ y/P{uJi)\ipe)- (36) 

Therefore only half the number of operators are required 
for the estimation experiment than when density matri¬ 
ces are used. 


VI. CONCLUSION 


V. DISCUSSION 

For our estimation scheme we utilized sequential mea¬ 
surements on a single quantum system instead of contin¬ 
uous measurements. A sequential measurement can be 
compared to a continuous measurement with the same 
measurement strength [7]. The measurement strength 
7m is a specific ratio of the strength of the individual 
measurements Ap and the time r in between measure¬ 
ments (cp. Section |n|. For continuous measurements the 
time steps r are made infinitely small, so for the same 
measurement strength many more integrations per qubit 
cycle are required than time steps needed for sequential 
measurements. The estimation of the dynamical param¬ 
eters of a quantum system is thus computationally more 
efficient with sequential measurements. 

The technique of monitoring the state and frequency of 
the Rabi oscillations of a system via sequential measure¬ 
ments is characterized by two timescales the level 

resolution time which defines the timescale on which 
the state evolves due to the measurement sequence (cp. 
Sec. 0 and Tfi which is the Rabi period, also char¬ 
acterizes the information gain due to measurement. In 
the case where the measurements dominate the evolution 
(Tm ‘C Tfj) the oscillations are modified and slowed down 


We have shown that it is possible to achieve full pro¬ 
cess tomography as well as real-time state monitoring by 
upgrading a known state estimation protocol. For this 
purpose we introduced one of the first applications of 
the hybrid quantum-classical formalism. Within this for¬ 
malism our estimated state is a hybrid state comprised 
of a probability distribution of the unknown parameters 
and a density matrix for the quantum state. Updating 
the quantum part of the hybrid state according to the 
measurement record induces an automatic update of the 
probability distribution according to Bayes’ Law. Nu¬ 
merical simulations were used to confirm the theory. 
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